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BEHAVIOR OF A SINGLE -TORN MAGNETIC COIL 


UNDER DYNAMIC LOAD CONDITIONS 

By M. D. Williams and G. K. Oertel 
Langley Research Center 


SUMMARY 


In a magnetic plasma compression experiment, a single-turn coil is sub- 
jected to magnetic pressure pulses. The energy developing the magnetic pressure 
pulses is stored in a 1-megajoule capacitor bank and delivered to the coil in a 
damped ringing electrical discharge. The kinetic pressures which can be devel- 
oped in the coil depend on the coil T s ability to contain these pressure pulses 
and that, in turn, is determined by the coil's mechanical design. Previously, 
all mechanical designs were based on calculations performed on the assumption 
that the peak current (and hence magnetic pressure) existed in the coil at all 
times during the discharge. This approach has a self-contained but unknown 
safety factor. In actuality, the coil undergoes a dynamic response throughout 
the discharge and is able to contain higher peak magnetic pressures than those 
predicted by static stress calculations because of an inertia term which is 
included in the dynamic response but is neglected in the static response. The 
dynamic response has been mathematically developed and the results have been 
compared with experimental data from a scaled-down model of the larger 
experiment . 


INTRODUCTION 


In many plasma experiments strong magnetic fields are generated for short 
times by the discharge of a capacitor bank through a single-turn coil. Often 
the discharge damps out long before the end of the first half-cycle of the 
mechanical oscillation of the coil. It was therefore thought that it might be 
quite possible for such coils to withstand much higher dynamic loads under 
these conditions than would be indicated by static stress calculations. 

Dynamic stress calculations were therefore carried out for a thin single- 
turn coil which is much longer than its radius. In addition to the length 
requirement, it was assumed that the amplitude of the mechanical oscillation is 
always much smaller than the radius of the coil and that the pressure amplitude 
is not affected by the mechanical oscillations. These assumptions are satis- 
fied in this experiment. The results apply to a range of mechanical oscilla- 
tion amplitudes with stresses ranging from zero to the elastic limit of the 



material. The elastic limit in turn determines what is herein defined as the 
threshold pressure. 

The predictions based upon the dynamic stress calculations were used in an 
experiment under conditions which are similar to those in a 1-megajoule magnetic 
compression experiment which is currently being designed and built. 

The results could have consequences for the maximum attainable magnetic 
pressures in coils and, according to the relation (see ref. 1, p. 53) 


p.. = § — = NkT = p, . 

M 2 kin 


it may then become possible to obtain higher temperatures and densities in the 
laboratory than were previously expected. (The formula given herein is the mks 
equivalent of the formula given in the cgs system in the reference.) This 
higher pressure capability could be important for the production of certain 
spectra of very highly ionized elements (mostly in the soft X-ray region) under 
controlled corona-like conditions as well as for the production of thermonuclear 
laboratory plasmas. 


SYMBOLS 


The rationalized mksAK system of units is used herein. 
Ai>Bi,Ci constants used in general solution 

B magnetic flux density 

E modulus of elasticity 

k Boltzmann constant 

l length of coil 

N particle density 

p pressure 

r(o) original inside radius of coil 

S yield strength 

T temperature 

t time 

u displacement 


2 



a 

damping constant 



constants 


7 

coupling constant 


S 

skin depth 


e 

strain 


n 

constant 


CD 

o 

V* 

CD 

H 

phase constants 


V 

permeability of coil 

material 

Vo 

vacuum permeability 


P 

mass density 


a 

conductivity of coil 

material 

T 

thickness of coil 


C0 o 

natural angular frequency of coil 

a>i 

angular frequency of 

electrical discharge 

Subscripts: 



dyn 

dynamic threshold 


e 

elastic 


kin 

kinetic 


M 

magnetic 


max 

maximum 


s 

static threshold 



Dots over symbols indicate derivatives with respect to time. 
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EQUATION OF MOTION 


The coil is considered to he a cylindrical ring with inside radius r( 0 ), 
thickness r, and length l . If p is the mass density of the material, E 
the modulus of elasticity, and u the displacement defined hy 

u = r(t) - r(0) 

in terms of the values of r at times t and 0, the equation of motion can 
he derived from a Lagrangian with the following terms (ref. 2 ): 


kinetic energy 
stress potential 

driving potential 


rtpTlr(0)&2 


-rtE-rl 


u c 


r(0) 
2 rtpZr( 0 )u 


where p is an internal driving pressure. The coil radius as a function of 
time r(t) may he replaced hy the coil inside radius r(0) in the driving 
term if the condition 


u « r(0) 


(la) 


or 


P « 


Et 

c( 0 ) 


(lh) 


is fulfilled. It is convenient to express the pressure as p = p ma> r f (t) , where 
f(t) is dimensionless and p mflY is independent of time. If, finally, the 
mechanical frequency o> 0 and the coupling constant 7 (an acceleration) are 
introduced hy 


and 



E 

pr 2 (0) 


_ p max 

' 2p T 


( 2 ) 


respectively, the equation of motion becomes 


u + o^^u = 2/f(t) 


( 3 ) 


k 



I ■ II ■nil llll I I II ■■III! IB ■■ 


If, in particular, the press-are is generated by a damped (damping con- 
stant a/2) oscillatory discharge with frequency o>j_ through the coil, then 
f(t) takes the form 

f(t) = e'^sin^GDjt (4) 

Also, if the condition 


per » 


1 


CD 2j2 


is fulfilled, that is, if the skin depth 5 is much smaller than the coil 
thickness t (p and a being the permeability and conductivity of the coil 
material, respectively), then 


P M 2p 0 

If 5 -> t, part of the magnetic flux could penetrate the coil material and var- 
iations of p^ would result along the coil. The additional requirement 

Z » 2r(0) could be eliminated in the experiment. 


EXPLICIT SOLUTION FOR SINGLE-TURN COILS 


Equation ( 3 ) in conjunction with equation (4) can be solved by the Laplace 
transform method (ref. 3)* With the definition 

P 2 - a 2 + cd 0 2 + 4-cd-j 2 

and the initial conditions u(0) = u(0) = 0, the following solution is obtained: 


u(t) = o 5- ie“ at - cos 0D o t + 7^- sin oo 0 t) 

+ cd qi 1 o J 


P^ - 1 6 cD Q 2 CDp 2 


,-at 


4clcd 2_ sin 2 co 2 _t - ^p 2 - ^l^) cos 2cD l^j 


+ 




COS CD 0 t 


a 

CDq 


P 2 sin 



(5a) 
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By combining terms, equation (5a) may "be "brought into the form 


where 


u(t) 




A-, = 1 

a 2 + oJq 2 


(5b) 



- I^q 2 ^ 2 

p 2 2 = 8cD! 2 - P 2 



U 1 


(t) = 


U>n 


: (“o 2 - Vl 2 ) 


(<*b 2 - V 2 )< 


-at 


CCt> 2 e -at cos 


2a>j_t + 4 <Dq_^cos a) 0 t 


(6) 


The threshold pressure can now be defined as that pressure which is just 
large enough to cause a permanent coil deformation . This threshold pressure is 
achieved if the elastic limit S is reached in the coil material (ref. 4). 

The static threshold pressure p s is related to the elastic limit by 



The strain e e which corresponds to the elastic limit is related to the thresh- 
old pressure by 


u e ^ Ps r (°) = e 
r(0) rE E 


(7) 


(which can he derived from equation ( 3 ) under static conditions). Assume that 
under dynamic conditions also the coil deforms permanently when the strain 
exceeds e e ; then in the limit of equation (6) the dynamic threshold pressure 
is given hy 


Pdyn - r(Q) 



«b 2 + oat_ 

Ir ^ 2 2 °* 


2 ( 1^1 


» 




» 


a) 


(8) 


This relation can he verified hy expanding equation (5a) in powers of a. If 
only first-order terms are retained, one arrives at equation (6) with (l - at) 
in place of the exponentials. The maximum elongation occurs when o^t » it. 

When 2 o>l » o^, the term with cos 2o>|_t will peak at least once near the peak 
of the third term; thus, cos 2a)j_t can he set equal to 1 to find the overall 
maximum value of u. 

It is not surprising that for a = 0 (undamped pressure, varying with 
sin^o^t) the dynamic and static threshold pressures are nearly equal in this 
limit. The average pressure for a = 0 is p mx ^2, and the coil oscillates 

about the position corresponding to this pressure with an amplitude which cor- 
responds to Pmax/ 2 ? thus, it reaches a displacement corresponding to Pmax 

during every cycle. With mechanical damping, equation (8) should include an 
additional factor of approximately 1 + itTj/o^ where r\ « o^, 2o^ is a 
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damping constant which would have come from an additional term +2tju in equa- 
tion (k) . In practice this correction is very small for most metals when 

T) « 


EXPERIMENT 


A short coil was tested with a small capacitor hank under stress condi- 
tions that are similar to those which will he encountered in a planned 
1-megajoule magnetic compression experiment. The experimental setup is shown 
in figure 1. The capacitor hank is connected to the coil hy high-inductance 
busses and can he switched with a single triggered air gap. Most of the induct 
ance is in these connections; only a few percent of the total inductance is in 
the coil. The most important parameters of the discharge circuit and the coil 
are given in table I; values obtained for the simulation experiment are com- 
pared, if possible, with those for the planned magnetic compression experiment. 

In scaling down to lower energies, it was important to keep the coil 
radius the same; therefore, the length had to be reduced considerably. Ordi- 
narily, this length reduction would mean a sizable increase in inductance, 
accompanied by an undesirable nonuniformity of the magnetic field. Both prob- 
lems could be solved by inserting an aluminum slug into the short coil with a 
gap between the slug and the coil of only about 3 percent of the coil radius. 
The skin effect excludes the high-frequency electromagnetic field from the slug 
this forces the field lines to be parallel to the coil and reduces the 



Figure 1.- Experimental setup. 
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TABUS I.- DISCHARGE CIRCUIT AND COIL PARAMETERS 


Capacitance, C, pF 

Voltage rating, V, kv* 

Angular frequency, o>]_, l/sec 

e-fold time, 2/a, sec 

Coil radius, r(0), m 

Coil length, Z, m 

Coil thickness, t, m 

Coil material 

Modulus of elasticity, E, psi 

Modulus of elasticity, E, N/m 2 

Yield strength, S, psi 

Yield strength, S, N/m 2 

Density, p, kg/m3 

Angular mechanical frequency, o) 0 , l/sec . . . . 


Simulation experiment 

22 

30 

1.38 x 105 

2.2 X 10“4 

6.2 X 10- 2 
1.46 X 10- 2 

1.6 X 10-3 

Electrolytic copper 

1.7 X 10T 

1.2 x 10 11 

2.9 x icA 
2 x 10 8 

8.9 x 103 
6.1 x lcA 


Magnetic compression 
experiment 

3000 

20 

1.25 X 105 


6.7 x 10 -2 
1.2 
3 X 10" 2 
Beryllium copper 

1.9 X 107 

1.3 x ion 

lO^ 

6.9 X 108 

9 x 103 

6.6 x 10^ 


inductance considerably. The magnetic pressure in the gap could be controlled 
by varying the charging voltage of the capacitor bank. 

In the first series of tests the general reaction of the coil to the pres- 
sure pulses was studied in an effort to determine the pressure range in which 
deformations could be expected. For the sake of safety, the first tests were 
made at bank voltages of 2.6, 3-5; and 4kV; in other tests the voltages were 
increased in increments of lkV. (The full pressure as in the planned experi- 
ment would be reached if the initial bank voltage were 55kV.) No changes were 
detected up to a bank voltage of l 8 kV when a small permanent deformation was 
first recorded. This bank voltage corresponds to a peak magnetic pressure of 
8.3 X 10^ N/m 2 (1200 psi). A further increase in the permanent deformation was 
observed after tests at still higher pressures. 

After about three tests at a particular pressure, the coil ceased to 
expand farther because of the increasing nonuniformity of the magnetic field 
in the gap, accompanied by a slight drop in pressure, and because of possible 
work-hardening effects in the coil material. 

Ordinarily, work hardening does not take place unless the elastic limit 
has been exceeded, which in turn would imply some permanent deformation. In 
order to reduce the possible effects of work hardening prior to the first 
observed permanent deformation, a second series of tests was made with a new 
coil, starting at a bank voltage of 9 ^V and increasing the voltage in incre- 
ments of 2 kV. The first permanent deformation was detected at a bank voltage 
of 17kV, corresponding to a peak magnetic pressure of J .6 X 1C)6 N/m 2 (1100 psi). 
Averaged over various diameters, a permanent change of 1.5 X 10-4 m (0.006 in.) 
in diameter was measured. Again, further permanent deformations occurred in 
tests at higher peak pressures. 
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The displacement u of the coil oscillating according to equation (5) for 
values of up and in this experiment is shown in figure 2 . The pressure 

corresponds to a bank voltage of 17 kV. For comparison, the strain computed 
from the approximate equation (6) is also given in the figure. If the same 

peak magnetic pressure ( 7-6 X 10 ^ N/m^) had been applied statically, the stress 
in the coil would have been 

S = P<3 = 7.6 X 10 6 X 10 ~ 2 = 285 X 106 N/m 2 

S T 

1.6 X 10 J 

From the static stress -strain characteristic of the coil material given in fig- 
ure 3 it can be seen that this value is much greater than the breaking point of 
the material. Alternatively, it can be shown that the permanent strain 
(0.0012) sustained by the coil dynamically corresponds to a peak pressure of 
approximately l 60 X 10^ N/m^ on the static curve. (The location of peak pres- 
sure and strain is at the intersection of the characteristic curve and a line 
parallel to its linear portion at a distance on the abscissa equal to the per- 
manent strain.) This fact indicates that the inertial term has accounted for 

* . , 285 - 160 

approximately 44 percent of the peak dynamic stress. 



Figure 2.- Displacement as a function 
of time. 



In/Sn. 


Figure 3*~ Static stress-strain 
characteristic of coil 
material. 


^Copper is stress-rate sensitive due to crystalline effects also. However, 
the good agreement herein between experiment and theory when only the inertial 
effect is used indicates that the crystalline effect is small. 


10 




The peak strain predicted by equation (5a) can he obtained from figure 2 
by dividing the value of u(t) by the coil radius. The maximum strain pre- 
dicted by equation (5a) is 0.0027* The maximum strain observed is 0.0026. The 
maximum strain lies slightly in the nonlinear portion of the static character- 
istic of the material. (This is requisite for a permanent strain.) In the 
linearized theory presented herein, E should be calculated from a linear 
curve which best fits this nonlinear segment. (This fit insures that the same 
energy per unit volume will be delivered to the coil as would be delivered 
under the nonlinear segment of the characteristic.) The nonlinear character- 
istic (fig. 3) can be represented very accurately by 


p = 1.48 x 10 1:L e - 4.1 x 10 15 e 2 + 2.73 x I0 1 5e3 


(0 % e < 0.014) 


The desired relation is 


p = Ee 

When the area between the two curves is equated to zero, E is found to be 

0.86 X 10^ N/m^. This value of E is used to calculate the maximum strain 
from equations ( 5 a) and ( 6 ). 

The error in the electrical quantities is of the order of 5 percent, or 
about 10 percent for the magnetic pressure (because of the square dependence 
on the current). The errors in yield strength and modulus of elasticity are 
likely to be only 1 or a few percent each. Thus, the comparison of experiment 
and theory should show agreement within 12 percent or better. 

According to equation ( 5 a) the maximum strain during the test which first 
indicated that yield strength had been exceeded was within k percent of the 
maximum strain determined from the characteristic curve and permanent strain. 
The agreement is based on the "best fit" calculated value of E. (Operating 
slightly into the nonlinear region of the characteristic curve means E is a 
nonlinear function of e.) Use of the best fit value of E provides a means 
of comparing experiment and theory. Normally, E from the linear portion of 
the characteristic curve would be used. 


CONCLUDING REMARKS 


Very good agreement between experimental and theoretical results is found 
for this experimental arrangement. The applicability of the simple relations 
to more general cases, such as thick coils, very high pulse frequencies, and 
long durations of the pressure pulse, remains to be tested. Within the frame- 
work of this formalism, it appears that the high frequency of oscillatory pres- 
sure pulses does result in significantly higher load capabilities in this 
experiment. Further, a short duration of the pressure pulses (i.e., a large 
damping factor) may lead to significantly higher load capabilities. Mechanical 
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damping is negligible for this purpose for most metals unless the pulse lasts 
much longer than one period of the mechanical oscillation. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September l 4 , 1965. 
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